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ABSTRACT 



Context. Axions are the natural consequence of the introduction of the Peccei-Quinn symmetry to solve the strong CP problem. All 

the efforts to detect such elusive particles have failed up to now. Nevertheless, it has been recently shown that the luminosity function 

of white dwarfs is best fitted if axions with a mass of a few meV are included in the evolutionary calculations. 

Aims. Our aim is to show that variable white dwarfs can provide additional and independent evidence about the existence of axions. 

Methods. The evolution of a white dwarf is a slow cooling process that translates into a secular increase of the pulsation periods of 

some variable white dwarfs, the so-called DAV and DBV types. Since axions can freely escape from such stars, their existence would 

increase the cooling rate and, consequently, the rate of change of the periods as compared with the standard ones. 

Results. The present values of the rate of change of the pulsation period of Gl 17-B 15 A are compatible with the existence of axions 

with the masses suggested by the luminosity function of white dwarfs, in contrast with previous estimations. Furthermore, it is shown 

that if such axions indeed exist, the drift of the periods of pulsation of DBV stars would be noticeably perturbed. 
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1. Introduction 

One of the long-standing problems of the standard model of par- 
ticles is the question of why charge and parity can be violated 
in weak interactions and not in strong interactions. This is the 
so-called strong CP problem. One of the possible solutions con- 
sists in the introduction of a new symmetry, the Peccei-Quinn 
symmetry or PQ symmetry, which automatically accounts for 
the problem (Peccei & Quinn 1977a; 1977b). The spontaneous 
breaking of this symmetry gives rise to the existence of a new 
particle called axion, which can couple to photons, electrons and 
r* baryons with a strength that depends on their mass and on the 
. £h ' s P ec ifi c wa y m which the PQ symmetry is implemented. The two 
K£ , most popular implementations are the KVSZ model (Kim 1979; 
j_i ■ Shifman et al. 1980), where the axions couple with photons and 
03 ; hadrons, and the DFSZ model (Dine et al. 198 1 ; Zhitnisky 1980) 
where they also couple to electrons. 

From the first moment that the problem was formulated, a 
huge effort aimed to directly detect these elusive particles, or at 
least to constrain their properties using astrophysical and cos- 
mological arguments, was undertaken. Among the experiments 
looking for a direct detection we mention ADMX and CAST — 
see Steffen (2009) for a synthetic description. Excellent reviews 
about the properties of axions and the corresponding astrophys- 
ical and cosmological searches can be found in Kim & Carosi 
(2008), Raffelt (2007; 1996; 1990), Khlopov (1999) and Turner 
(1990). 



White dwarfs represent the final evolutionary stages of low- 
and intermediate-mass stars. Since they are strongly degener- 
ate and do not have nuclear sources, their evolution is just a 
gravothermal process of cooling. The simplicity of these objects, 
the fact that the physical processes necessary to understand them 
are well identified, although not always well known, and the im- 
pressive existing observational background make these stars ex- 
tremely useful laboratories to test new ideas in physics (Isern 
& Garcfa-Berro 2008). In particular, Isern et al. (2008; 2009) 
found that the best fit to their luminosity function is obtained 
when axions of the appropriate mass are included in the calcula- 
tions. We note that in this case only DFSZ axions were taken into 
account, since electron bremsstrahlung is the dominant process 
in white dwarf interiors. 

The axion emission rate under these conditions is given by 
(Nakagawa et al 1987; 1988): 



e ax = 1.08 x 10 23 ^p — T*F(T,p) erg/g/s, 

47T A 



(1) 
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where F takes into account the Coulomb plasma effects, Tj is the 
temperature in units of 10 7 K, Z and A are the atomic and mass 
numbers of the plasma components, respectively, and g aee is the 
strength of the axion-electron Yukawa coupling, which is not 
defined by the theory and determines the total axion luminosity 
of the white dwarf. The value of this parameter that best fits the 
luminosity function of white dwarfs is 1 . 1 x 10~ 13 , but variations 
of a factor two are still compatible with the observations (Isern 
et al. 2008; 2009). 
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Notice that strictly speaking, the shape of the white dwarf 
luminosity function can be equally well fitted invoking the exis- 
tence of any hypothetical light boson or axion-like particle able 
to couple to electrons with a strength g aee ~ 10~ 13 . However, 
since axions have a solid theoretical justification, we will use 
this term all along the paper. The relationship between this cou- 
pling constant and the mass of the axion is given by 



2.8 x 10 



_ 14 m ax cos 2 /? 
lmeV ' 



(2) 



where cos 2 /? is an essentially free parameter that probably takes 
a value ranging from 0.5 to 1 and comes from the ratio of 
the two Higgs fields that characterize the DFSZ axion models. 
Therefore, the axion mass necessary to fit the white dwarf lumi- 
nosity function mass is m ax cos 2 /? ~ 4 meV where, as previously 
mentioned, variations of ~ 2 meV are compatible with the obser- 
vational uncertainties. For simplicity we will assume cos 2 /? = 1 
and will express our results in terms of the mass of the axion, as 
it is usually done. 

Taken at face value, the result of Isern et al. (2009) can be 
considered as the first indirect proof that axions of the DFSZ 
type indeed exist. However, before accepting this asseveration 
some conditions have to be fullfiled: i) Different and indepen- 
dent evolutionary codes must provide the same result, ii) there is 
no conventional way to account for the presently observed slope 
of the white dwarf luminosity function, iii) there is no contradic- 
tion with other well established astrophysical or cosmological 
phenomena, and iv) there are some phenomena that are better 
explained if the axion hypothesis is adopted. Concerning this 
last point, it is important to realize that the value of the mass of 
the axion found in Isern et al. (2009) is close to the upper limit 
obtained from astrophysical constraints. This means that axions 
with this mass will not dramatically change our current under- 
standing of stellar evolution, but instead some subtle changes 
that accurate observations could detect and quantify would be 
apparent. The case of the white dwarf luminosity function is 
just an example. The drift of the pulsation periods of variable 
white dwarfs, the topic of this paper, is another one that will be 
discussed in the following sections. We note however that this 
is not the case of core collapse supernovae, where axions with 
this mass can drain an important amount of gravitational energy 
and modify the process of collapse/explosion. Unfortunately, the 
nucleon-axion interaction, which is dominant in supernovae, is 
not well understood and deserves additional work that is out 
of the scope of this paper. Finally note as well that if cos 2 /? 
was small, the mass of the axion would be higher than the val- 
ues quoted here and, consequently, the interaction axion-photon 
would be strengthened. This could be used to rule out our results 
or to constrain the values that cos 2 /? can adopt. 



2. Results and discussion 

During the process of cooling, white dwarfs experience several 
phases of pulsational instability powered by the /(-mechanism 
and the "convective driving" mechanism (Winget & Kepler 
2008). Depending on their location in the Hertzsprung-Russell 
diagram and their atmospheric composition they are called pul- 
sating PG1159 or GW Vir stars, which includes stars with a 
surrounding nebula (the PNNV stars) and stars without nebula 
(the DOV stars) at r eff ~ 70,000 - 170,000 K and He-, C-, 
and O-dominated atmospheres, DBV or V777 Her stars, with 
r e ff ~ 25,000 K and He-dominated atmospheres, and DAV or 



Fig. 1. Evolution of the measurement of the period drift of 
Gl 17-B15A. Lines represent the theoretical values obtained by 
Corsico et al. (2001) — solid line — and Bischoff-Kim et al. 
(2008) — dashed-dotted line (for the case of a thin envelope) and 
dashed line (for the case of a thick envelope). The nominal un- 
certainties of these models are 1.15, 0.26 and 0.17 respectively, 
in the units of the figure. 



ZZ Ceti stars, with T e tf ~ 12,000 K and H-dominated atmo- 
spheres. All of them display multiperiodic pulsations with pe- 
riods in the range of 100 to 1000 s, although PNNV stars can 
pulsate with periods as long as 3000 s. The length of these peri- 
ods clearly indicates that these objects are experiencing g-mode 
non-radial pulsations, where the main restoring force is gravity 
— see Winget & Kepler (2008) for a recent review. 

One of the main characteristics of g-mode pulsations is that 
they experience a secular drift that can be used to test the white 
dwarf cooling theory. This secular drift can be approximately 
described by 



c/lnn 

dt 



dlnT d\nR 
-a — - — + b- 



dt 



dt 



(3) 



where a and b are constants of the order of unity that depend on 
the details of the model, and R and T are the stellar radius and 
the temperature at the region of period formation, respectively 
(Baglin & Heyvaerts 1969; Winget et al. 1983). This equation 
reflects the fact that as the star cools down, the degeneracy of 
the plasma increases, the buoyancy decreases, the Brunt-Vaisala 
frequency becomes lower and, as a consequence, the spectrum 
of pulsations gradually shifts to lower frequencies. At the same 
time, since the star contracts, the radius decreases and the fre- 
quency tends to increase. In general, DAV and DBV stars are 
already so cool that the radial term is negligible and the change 
of the period of pulsation can be directly related to the change 
in the core temperature of the star. The timescales involved are 
of the order of ~ 10" 11 s/s for DOVs, ~ 10~ 13 to ~ 10~ 14 s/ for 
DBVs and ~ 10~ 15 to ~ 10~ 16 s/s for DAVs. The measurement 
of such drifts is a difficult task but, as was proved, it is feasible 
(Winget & Kepler, 2008). 

These properties allow us to build a simple relationship 
(Isern et al. 1992; Isern & Garcfa-Berro 2008) to estimate the 
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Table 1. Axion luminosities and axion masses necessary to ac- 
count for the observed drift of the period of pulsation of Gl 17- 
B15A. 



Model 


L ax /10 ju (erg/s) 


ra ax cos 2 /? (meV) 


C 


0.6(2.8/0.0) 


2.8 (6.0/0.0) 


BK (thick) 


6.8(8.7/4.9) 


9.3 (10.6/7.9) 


BK (thin) 


2.2(3.5/1.0) 


5.5 ( 6.9/3.7) 



influence of an extra sink of energy, like axions, on the period 
drift of variable white dwarfs: 

-^model n m0( ] e i Tmodel 

where the suffix "model" refers to those models built using stan- 
dard physics. 

G117-B15A is a ZZ Ceti star (McGraw & Robinson 1976) 
with a dominant period of pulsation of ~ 215 s. This star, in 
an exceptional observational effort, has been observed for more 
than 30 years to obtain the secular evolution of the pulsation 
period. Kepler et al. (1991a; 1991b), after ~ 15 years of obser- 
vations, found a secular drift of n o b s = (12.0 + 3.5) x 10~ 15 s/s 
that was more than a factor of two larger than the expected value 
in the case of a normal carbon-oxygen white dwarf. In order to 
account for this discrepancy, Isern et al. (1992) proposed the in- 
troduction of an extra cooling term due to axions. The mass of 
the axions necessary to fit this discrepancy was 8.5 meV. Later 
on, the analysis of additional ten years of observations provided 
a value of fl obs = (2.3 + 1.4) x 10~ 15 s/s (Kepler et al. 2000), 
closer to the value predicted by the standard theory of evolution 
of white dwarfs. A detailed seismological analysis (Corsico et 
al. 2001) found a drift ri mode i = (3.9 + 2.3) x 10~ 15 s/s, where 
the main sources of uncertainty were due to the mode identifica- 
tion, the mass of the star, the chemical profile and the effective 
temperature. This value suggested that it was no longer neces- 
sary to invoke axions to account for the observed cooling rate of 
white dwarfs, but at the same time the uncertainties were large 
enough to prevent ruling out their existence. It is important to 
realize here that the values obtained with Eq. © agreed well 
with those obtained from the accurate numerical calculations of 
Corsico et al. (2001). A new analysis including additional five 
years of observations gave II = 3.57 + 0.82 x 10~ 15 s/s (Kepler 
et al. 2005), a value that overlaps with the theoretical values ob- 
tained by Corsico et al. (2001) and is marginally compatible with 
the axion luminosity obtained if the axion mass suggested by 
the luminosity function of white dwarfs, ~ 4 meV is adopted. 
Finally, in a recent release that includes yet another five years of 
observations, Kepler (2009) obtained fl = 4.77 +0.59 x 10 -15 s/s. 
Adopting the proper motion correction of Kepler et al. (2005), 
we obtain IT = 4.07 + 0.61 x 10~ 15 s/s, a value that is consis- 
tent with the existence of the extra cooling term suggested by 
the luminosity function. Figure [TJ displays the evolution of the 
measured value of II with time. The jumps between the differ- 
ent observational values reflect the introduction of corrections, 
like the proper motion correction, and the fact that the phase of 
the periodicity is affected by a jitter caused by low-amplitude 
modes, contamination by G117-B15B, the companion star, use 
of different telescopes and apertures. . . that are not well under- 
stood yet (Kepler, private communication). 

TableQJdisplays the axion luminosity and mass necessary to 
fit the observations when several white dwarf pulsational models 
are adopted (the values in brackets correspond to the observa- 
tional uncertainties). The values of the first row, labeled model 



C, were obtained using Eq. ©, the model of G117-B15A that 
best fits the data (Corsico et al. 2001) and the previously men- 
tioned axion emission rates of Nakagawa et al. (1987; 1988). It is 
important to realize here that in a completely independent anal- 
ysis, Bischoff-Kim et al. (2008) identified two possible astero- 
seismological models of Gl 17-B15A, one with a thin hydrogen 
envelope and the other with a relatively thick hydrogen enve- 
lope, with II = 2.98 + 0.17 x 10~ 15 s/s and 1.92 + 0.26 x 10~ 15 
s/s, respectively. Both values are smaller than the measured ones 
and indicate that an additional sink of energy is necessary. Table 
[TJ shows the luminosity and the mass of the axions necessary to 
account for the most recent data of Kepler (2009) assuming a 
temperature of the core of T c = 1.2 X 10 7 K in both cases. 

Recently, Mukadam et al. (2009) analyzed R 548 (ZZ Ceti 
itself) and found that its mass lies between 0.55 and 0.58 Mq, 
while the mass of the hydrogen envelope is log Mu — -5.0 and 
that the best prediction for the rate of change of the period is II = 
(8.3 + 0.046) x 10~ 15 s/s. These results are fully consistent with 
the values obtained for Gl 17-B15A and with those displayed in 
Table [TJ 

The pulsation properties of DBV stars can provide an addi- 
tional test to the axion hypothesis. Since these stars have a hotter 
core than that of DA white dwarfs, the neutrino luminosity can 
no longer be neglected, and the influence of axions on the pul- 
sation period drift is smaller than that expected for DA white 
dwarfs. For instance, a DB white dwarf of M ~ 0.59 Mq, an 
effective temperature T e ff ~ 25 200 K and a core temperature 
log T c = 7.61 has a luminosity of L ~ 6.3 x 10~ 2 Lq and a neu- 
trino luminosity L v ~ 5.1 x 10~ 2 Lq. If an axion mass of 5 meV is 
assumed, the axion luminosity amounts L ax = 6.0 x 10 2 Lq and 
the corresponding secular change of the pulsation period would 
be 



1A no-ax 

Note that because of the strong dependence of the different lumi- 
nosities on the core temperature, this prediction is very sensitive 
to the position in the instability strip and on the particular param- 
eters of the star under study. The only case of DB variable white 
dwarf that is currently being studied is EC20058-5234 (Sullivan 
2009), but unfortunately a rate of period change is not yet avail- 
able for this star. 

3. Conclusions 

It has been shown that white dwarf variable G117-B15As 
present value of the secular rate of change of the period of pul- 
sation, 4.07 + 0.61 x 10~ 15 s/s (Kepler 2009), is consistent with 
the predictions of the theoretical models, as was anticipated by 
Isern et al. (1992), if an additional source of cooling like ax- 
ion emission is included. This result is corroborated by the com- 
pletely independent analysis of the same star done by Bischoff- 
Kim et al. (2008), using the measured rate of period change 
(3.57 + 0.82 x 10 -15 s/s) obtained previously (Kepler et al. 2005). 
This means that the conclusion of Corsico et al. (2001) that it 
was unnecessary to introduce an additional cooling source is no 
longer valid when the new observational data are taken into ac- 
count. This is an important result, as asteroseismological obser- 
vations of white dwarfs seem to give additional and independent 
support to the claim of Isern et al. (2008; 2009) that the white 
dwarf luminosity function is better fitted if an additional cooling 
provided by axions with a mass of a few meV is included in the 
calculations. 
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It has also been shown that since the core temperatures of 
DBV stars are higher, the secular drift of the pulsation period of 
these stars is strongly enhanced by the axion emission. In par- 
ticular, if axion emmisivity is computed using the mass obtained 
from the best fit to the luminosity function of white dwarfs and 
included in the calculations, the secular rate of period change 
increases by ~ 50%, although the exact value depends on the 
core temperature and, consequently, on their position in the DB 
variability strip. 

Finally, it is worthwhile to emphasize the importance of con- 
tinued obervations of G117-B15A and R 548. It is also impor- 
tant to enlarge the number of stars with observed period drifts 
and the measurement of these drifts in hotter variables like DBV 
white dwarfs. Finally, it is also important to improve the astero- 
seismological models. All these actions would undoubtely con- 
tribute to solve the problem of axions. 
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